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We present a new method for the solution of the Schrodinger equation applicable to problems of 
non-perturbative nature. The method works by identifying three different scales in the problem, 
which then are treated independently: An asymptotic scale, which depends uniquely on the form of 
the potential at large distances; an intermediate scale, still characterized by an exponential decay of 
the wave function and, finally, a short distance scale, in which the wave function is sizable. The key 
feature of our method is the introduction of an arbitrary parameter in the last two scales, which is 
then used to optimize a perturbative expansion in a suitable parameter. We apply the method to 
the quantum anharmonic oscillator and find excellent results. 
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In this Letter we present a new method to find ap- 
proximate solutions to eigenvalue problems in Quantum 
Mechanics. In particular, this method is applicable to 
problems where ordinary perturbation theory generates 
divergent asymptotic series. Actually, improving the con- 
vergence of the standard Rayleigh-Schrodinger perturba- 
tive expansion has been the subject of many studies in the 
past (see, e. g., Refs. an d references therein) 

and many variants of "optimized expansions" have been 
proposed. 

Among them we would like to single out the Linear 
Delta Expansion (LDE) ||, which is at the base of our 
developments. The LDE has been extensively applied in 
many different settings with varying degrees of success. 
For example, in Ref. jQ| it has been used to analyze dis- 
ordered systems. In Ref. Q it has been applied to study 
the slow roll potential in inflationary models. Pinto and 
collaborators have applied it to the Bose-Einstein con- 
densation problem |f, tnc 0{N){<j> 2 )\ d model 0, the 
Walecka model ^3 an d to the 4 theory at high tem- 
perature More recently the LDE has been applied 
with success to the study of classical nonlinear systems 
[l2L . Detailed references on the method can be found 
in these works. 

Here, we illustrate the method we propose by applying 
it to the quantum anharmonic oscillator, which is the 
usual benchmark to test any non-perturbative method. 

The cornerstone of the method is the identification of 
three different scales in the problem, which give rise to 
different behaviors of the wave function. Keeping in mind 
the standard separation of the Hamiltonian in a unper- 
turbed piece (the harmonic oscillator) and a perturbative 
one (the anharmonic term), one can recognize that at 
very large distances the wave function assumes its asymp- 
totic behavior: This is completely determined by the an- 



excited) states; at intermediate distances the wave func- 
tion still decays exponentially, but now governed by the 
harmonic term; finally, there is a short distance scale, in 
which the wave function is sizable. 

The key feature of our method is the introduction of 
an arbitrary parameter in the last two scales, which is 
then used, in the LDE spirit, to optimize a perturbative 
expansion in a suitable parameter. 

Consider the Schrodinger equation 
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%pn{x) = E n 1p n (x), (1) 



where /i is the anharmonic coupling. The asymptotic 
behavior of tp n (x) in the region of large x is determined 
by substituting the ansatz i[) n {x) oc e _7 ' 3: ' P into Eq. 0). 
One obtains p = 3 and 7 = \J fxm/2/2>'h. 

In order to make the three scales explicit in the wave 
function, we write 
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where the exponential takes care of the correct behavior 
in the limit \x\ — > 00. Notice that the quadratic term 
in the exponential docs not affect the behavior at large 
distances, but dominates at scales where |x| < f3/j. Here 
/3 = m\Juj 2 + £l 2 /2h is the coefficient of a harmonic os- 
cillator of frequency (l = \/ lo 2 + £l 2 , where f2 is an arbi- 
trary parameter introduced by hand; j£ n is a well-behaved 
function, which fulfills the equation 
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Equation J5J) has been introduced in order to single out 
three different regimes in the wave function: The purely 
asymptotic regime {\x\ — > oo), where the cubic term 
in the exponential dominates; the intermediate regime, 
where \x\ is not yet asymptotic but large enough to ex- 
pect the wave function to be exponentially damped; the 
regime of small \x\ where the physics is all contained in 
the £'s. The last two regimes will display a dependence 
upon the arbitrary frequency Q, although in a quite dif- 
ferent fashion (in fact, the intermediate regime displays 
a truly non-perturbative dependence upon Q). In the 
limits (/i, O) — one obtains the equation for the har- 
monic oscillator of frequency lu, which admits polynomial 
solutions (the Hermite polynomials). 

It is worth stressing that the energy E n in Eq. Q) is 
still the true energy, since no approximation has been 
used to derive this equation. 

We observe in Eq. J2J) that both the wave function 
£n( x ) an d the energy E n depend in some nontrivial way 
upon the anharmonic coefficient (i. On the other hand, 
the dependence of ip n (x) and E n upon the arbitrary fre- 
quency f2 is only fictitious, since this parameter does not 
appear in the original equation |T]). Nonetheless, we will 
show that Q can be used to generate an efficient expan- 
sion for the solution of Eq. . 

Indeed we rewrite Eq. @ as 



minimized byapplying the Principle of Minimal Sensitiv- 
ity (PMS) llj, i.e. by requiring that a given observable 
O (the energy, for example) be locally independent of fl: 
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80. 



0. 



(6) 



We will illustrate the method by explicitly showing the 
first two orders, although we have obtained results up to 
the eighth order. To lowest order Eq. reduces to the 
equation of a harmonic oscillator of frequency Q, whose 
solutions are the Hermite polynomials, 



while the energy eigenvalues are given by 

r 



E n Q 



HQ n 



n = 0, 1,... 
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To first order we have the equation 
2m£lx 
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where the left hand side of Eq. corresponds to the 
equation for a harmonic oscillator of frequency Q. Fol- 
lowing the spirit of the LDE, we have introduced a pa- 
rameter 8, which is going to be used as a power-counting 
device: When 8=1, Eq. reduces exactly to Eq. ©. 

Although 8 is not a small parameter wc will treat the 
right hand side of Eq. Q as a perturbation, writing down 
the following expansions: 



Ux) = Y, 5 ^nj{x), E n =J2 S j E, 
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Combining Eqs. 10} and JSJ, one can generate a hierarchy 
of equations, corresponding to the different orders in 8. 

Since we are doing perturbation theory, all the results, 
to any finite order in the expansion, will depend upon the 



Although such equation is valid for any state of the AHO, 
for illustrative purposes we will now consider only the 
ground state, for which ^oo( a; ) = 1- Then, for the case 
n = 0, the solution of Eq. © is a polynomial of order 3 
and can be cast in terms of unknown coefficients as: 



£oi(aO 



ao + a\x + a-ix 2 + a^x 3 . 



(10) 



The coefficient ao is not determined by Eq. 10 and we 
impose ao = 0, since it corresponds to the same func- 
tional form of £oo 0- By substituting this polynomial 
into Eq. © one gets the coefficients 



a>l = 0; a 2 
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and the energy 
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Therefore, up to first order, we get the energy 
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It is important to notice that the wave function obtained 
up to first order does not have nodes, a desirable result 
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The PMS to this order yields the solution f2 
the corresponding energy 



E 
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Here we would like to stress a point: We could have 
obtained a solution similar to that of Eq. (|1(J|) by di- 
rect application of the Raylcigh-Schrodingcr perturba- 
tion theory to the wave functions of the harmonic oscil- 
lator. However, had we applied the LDE directly to the 
Raylcigh-Schrodingcr expansion, we would not have been 
able to reproduce the correct asymptotic behavior of the 
wave function Q. 

The extension of the method to include higher orders 
is straightforward and the details will be presented else- 
where [l5j . Here we just present the expression for the 
energy up to third order 
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The optimal value of f2 is then obtained by using the 
PMS. The expression for fl in the general case is lengthy 
(see 0) and we only write here its asymptotic value in 
the limit of large [i (from now on we assume Ti = m = 
lu = 1): 
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(16) 



We can then substitute this expression in Eq. (|15|) and ex- 
tract the asymptotic behavior of the energy of the ground 
state to third order: 
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The asymptotic behavior of the ground state energy of 
the AHO has been studied in [l6j, by using the exact 
Raylcigh-Schrodingcr perturbation coefficients as an in- 
put. We use this calculation as a reference to compare to 
our results. Following the notation set in |l6j we write 
the energy as 
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and using this formula in Eq. I|17|) we obtain the coeffi- 
cient «o to third order in our expansion to be 



(3) 
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0.627, 



(19) 



which falls within 6 % of the true value. We have also 
computed ao up to eighth order, where the agreement is 
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FIG. 1: The coefficient ao in the asymptotic expansion of 
the ground state energy (see text) at different perturbative 
orders, up to order 8. The solid line is the exact result [l6|. 




FIG. 2: Ground state energy of the AHO as a function of 
to different perturbative orders, assuming h = m = u = 1. 



In Fig. |21 we plot the energy of the AHO as a function 
of [i. The solid bold curve is the exact result obtained 
by numerically solving the Schrodingcr equation for the 
AHO. The other curves are the energies computed at dif- 
ferent perturbative orders, up to fifth order (from bottom 
to top). We notice from the plot that the perturbative 
result is always below the true result; this is to be con- 
trasted with a variational calculation, where the exact 
result would be approached from above. We do not know 
whether this is a general property. 

In Fig. |21 we show the energy of the AHO as a function 
of the parameter f2, for fi = 5. The different curves cor- 
respond to different orders in perturbation theory (from 
second up to fifth order). We observe that the energy 
at different orders in the perturbative expansion devel- 
ops local minima and maxima. By applying the PMS we 
select the extremum in which the energy is flatter, i. e. 
where the dependence upon f2 is smaller. 
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displayed results up to a few orders in the perturbative 
expansion, it is very easy to push the calculation to any 
order, since the method only requires the solution of al- 
gebraic equations order by order. 

We are currently working on the application of our 
method to the calculation of the excited states of the 
AHO, to more general anharmonic potentials and to the 
double-well potential. 

P. A. and A. A. acknowledge support for this work to 
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also acknowledges Conacyt grant no. C01-40633/A-1. 



FIG. 3: Energy of the AHO as a function of SI, at different 
perturbative orders (% = m = U) = 1) and keeping fj, fixed 
(/* = 5). 




FIG. 4: Wave function for the ground state of the anharmonic 
oscillator, assuming m = h = u = 1 and fi = 200. The solid 
line is the exact result whereas the dashed (dotted) line is the 
result obtained with our method to third (fourth) order. 



well for the wave function. Indeed, in Fig. 0] we plot the 
wave function of the ground state of the AHO. The solid 
line is the exact (numerical) result, whereas the dashed 
and dotted lines refer to the wave function obtained by 
applying our method to third and fourth order respec- 
tively. We see that it works extremely well, even for large 
values of [i (/i = 200). Although we have not imposed 
that the wave function be flat at the origin, we observe 
that the method naturally provides this result. 

Finally, we would like to stress that, although we have 
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We are considering only the region x > 0. The other 
region will be obtained by using the symmetry properties 
of the wave function. 

This is somewhat analogous of the procedure employed 
in the Lindstedt-Poincare method to get rid of "secular 
terms" in the solutions fl2Tll3l . 



